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SUMMARY 


Closed lunar orbits are envisaged in lunar mission programs. The 
study described herein was undertaken to obtain an appreciation of the 
relevant fuel consumption requirements. The retrograde impulses nec- 
essary for establishing the orbits were assumed to occur at the point 
of closest approach of the main earth-moon trajectory; this point, 
designated as the arrival position, was restricted to a lunar altitude 
of ^>,000 nautical miles or less. The orientation of the arrival posi- 
tion vector relevant to any coplanar radius vector is not constrained, 
however, and similarly the scalar value of the arrival velocity is 
unrestrained. 

Since the arrival altitude is restricted to 5 >000 nautical miles 
or less, the perturbing accelerations of the earth and sun are suffi- 
ciently small that the vehicle and moon essentially comprise an isolated 
two-body system; this is discussed in the report. 

Retrograde velocities are determined for any required pericynthion 
position. If the pericynthion orientation requirement is relaxed then 
a smaller retrograde velocity is in some cases possible. A comparison 
between minimum retrograde velocities and retrograde velocities neces- 
sary for stipulated pericynthion positions is given. Arrival, veloci- 
ties are correlated with feasible earth departure conditions. 

The equations developed for determining retrograde velocities for 
desired pericynthion positions are considered useful for estimating 
essential data for the preliminary planning of lunar missions. Some 
graphical representation is included herein for immediate familiariza- 
tion with possible conditions. 


INTRODUCTION 


The study described in this report was initiated by the desire for 
a preliminary assessment of various problems associated with local, lunar 
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orbits as envisaged in lunar mission programs. Since it was known that 
the lunar gravitational field would be dominant during low-altitude 
orbits and that two-body solutions would be applicable, it was decided 
to use a closed-solution approach. When satellites of a dominating 
gravitational field are studied, velocities relative to the nonrotating 
coordinate system, with origin at the gravity source, are treated as 
inertial in order that Newtonian laws may apply. As an example, planet 
velocities relative to the sun are considered inertial and yet the sun 
is thought to be moving in space. For a lunar satellite, a perfectly 
elliptical orbit cannot be achieved due to the movements of the sun, 
earth, and moon relative to the vehicle; however, for low-altitude lunar S 

orbits, the lunar gravitational field is dominant; hence, near elliptical 1 

orbits may occur. Stability checks of the ellipticity were examined 0 

from the output of a digital program employing the Encke technique of 
integration with origin at the moon center. This was considered to be 
sufficiently accurate and the checks confirmed the validity of the 
approach taken. 

Contained within this report is a method for determining the retro- 
grade velocity necessary for establishing orbits with specific charac- 
teristics. It is anticipated that the plane of the vehicle f s lunar 
arrival velocity will, by guidance impulses, nearly coincide with the 
plane of the desired resultant orbit; hence, one constraint of this 
study was that both the arrival velocity and resultant orbit be coplanar. 

This report also indicates the correlation of arrival velocity with 
specified earth departure conditions by means of a restricted three- 
body mathematical model. 

It is possible in many cases to reduce Ihe magnitude of the retro- 
grade impulse if the restrictions on the perj cynthion orientation are 
relaxed. Therefore, a method is presented whereby the minimum retro- 
grade impulse can be determined for an orbit with a specified peri- 
cynthion radius with no constraint on orientation. The methods of 
determining minimum retrograde velocities anc. retrograde velocities 
for required response are thought to be us eft 1 in that they point the 
way to a comprehensive quantitative survey program for assessing fuel 
requirements for lunar orbits. 


SYMBOLS 

a semimajor axis of orbit, ft 

C Jacobian constant, ft^/sec^ 


D 


distance from center of earth tc center of moon, ft 



Naperian logarithm base 

universal gravitational constant, ft^-lb/sec^ - 

gravitational acceleration at earth's surface, ft/sec 2 

altitude, nautical miles 

altitude at apocynthion, nautical miles 

altitude at pericynthion, nautical miles 


specific impulse of fuel, sec 

mass of earth plus mass of moon, slugs 

mass of earth, slugs 

mass of* moon, slugs 

mass of* lunar vehicle, slugs 

distance from moon center to vehicle position, ft 
apocynthion distance from moon center, ft 

distance from barycenter to vehicle position, ft 

distance from earth center to vehicle position, ft 

pericynthion distance from moon center, ft 


distance from earth center to barycenter, ft 
distance from moon center to barycenter, ft 
time, sec 


velocity referred to the rotating coordinate system, 
the barycenter, ft/sec 


origin 


arrival velocity referred to the nonrotating 
system, origin the moon center, ft/sec 


coordinate 
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V e 

Vi 

V m 

Vo 

V oa 


velocity referred to the nonrotating coordinate system, 
origin the earth center, ft/sec 

velocity referred to the nonrotating coordinate system, 
origin the barycenter, ft/sec 

velocity referred to the nonrotating coordinate system, 
origin the moon center, ft/sec 

orbit velocity referred to the nonrotating coordinate 
system, origin the moon center, ft/sec 

orbit apocynthion velocity referred to the nonrotating 
coordinate system, origin the moon center, ft/sec 


Vop 


V 


R 


orbit pericynthion velocity referred to the nonrotating 
coordinate system, origin the noon center, ft/sec 

retrograde velocity, ft/sec 


X,Y,Z rotating coordinates, origin the barycenter 

X e ,Y e ,Z e nonrotating coordinates, origin the earth center 

x i> Y i> z i nonrotating coordinates, origin the barycenter 


x mJ Y m> z m 


7 

8 
6 


€ 

9 




M 


nonrotating coordinates, origin :he moon center, X m and 
Ym in the plane of the vehicle's orbit 

angle between vehicle velocity vector and local horizontal, 
deg 

ratio of fuel mass to gross mass at commencement of 
thrusting for entry into local lunar orbit 

overall ratio of fuel mass to gross mass required for both 
entry into and exit from the local lunar orbit at the 
same point of orbit 

eccentricity 

angle between r and rp at instant of retroimpulse 
burnout , deg 


S 

1 

0 
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s 

1 

0 



£ included angle between and Vr, deg 

<jd rotational velocity of earth-moon system, radians/sec 

A dot over a symbol indicates differentiation with respect to time. 


DISCUSSION 


Arrival Velocity 

In lunar mission programs, it is assumed that one of the mission's 
main objectives will be to make a close survey of the moon's terrain. 

For a detailed survey, this requires the vehicle to establish an orbit 
about the moon which will necessitate the application of a retrograde 
impulse. It is assumed for this study that the impulse will be applied 
at the instant of closest approach to the moon when the vehicle's velocity 
vector is normal to the extension of the moon's radius. A conception of 
the complete mission is shown in figure 1. 

The vehicle's lunar arrival velocity and position can be correlated 
with the earth insertion velocity and position by reference to a restric- 
ted three-body mathematical model. In this model the moon and earth are 
considered to rotate with constant radii and constant angular velocity 
about their common center of mass. The earth and moon are considered 
as point masses and the vehicle's mass is regarded as infinitely small 
in comparison. Further information concerning the characteristics of 
the restricted three-body mathematical model is found in reference 1. 

The restricted three-body equations of motion as given in the 
rotating coordinate system with the barycenter as origin are 


x . * 2ffi Y - 01(1 - "K* - x i) . °*( x - *?) 

r~„3 r 3 


Y = o^Y - 2 cjdX 


GM(1 - q)Y 



GMgY 

r3 


( 2 ) 



(3) 


-GM(l - n)Z GNLlZ 



where (x - X]_) is the distance along the X coordinate from the earth 
center to the particle and ^X - X 2 ) is the distance along this coordinate 
from the moon center to the particle. 

By writing 


W(X,Y,Z) = i o^(x 2 + y2) + (4) 

2 \ / r ep r 

then the equations of motion become 


X = — + 2ooY 

(5) 

8x 


Y = — - 2oiX 
dY 

(6) 

z - * 

(7) 


Sz 


By multiplying equations ( 5 ) to ( 7 ) by 2X, 2Y, and 2Z, respectively , 

adding, and integrating, Jacobi's integral :.s obtained as 


• P • 2 • 2 

A + Y + Z 


= 0 ? (x 2 + Y 2 } 


2GM(: 


jO + 2GMp 


- C 


ep 


( 8 ) 


or 


C 


o 2 + 2(1 - 4)GM + i-^GM _ y2 
b P r 


(9) 


where r-^p is in the earth-moon plane. 
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If the earth departure velocity and position are known, the integra- 
tion constant C can he determined. Once C is determined, it is 
possible to calculate the scalar value of the lunar arrival velocity for 
any given arrival position. The transformation of the velocity in the 
rotating coordinate system, origin the barycenter, to the velocity in 
the nonrotating coordinate system, origin the moon center, is 



1 The relationship between the velocities in the different coordinate 

3 systems is illustrated in figure 2. A comparison between earth departure 

velocities and lunar arrival velocities is shown in figure 3. The cor- 
relation between arrival velocity and earth departure velocity is shown 
for the condition of minimum earth-moon distance. The vehicle's arrival 
velocity is in the ^-Zo x direction. 

In acquiring a quantitative feel for lunar arrival velocities as 
a function of earth departure velocities for a given earth departure 
position, it is natural to think in terms of V e rather than V. A 
knowledge of the value of V does not immediately give the value of V e - 
however, for the departure position given in figure 3, it is known that 6 * 
y will be in the range of -2° to 25 ° in the case of coplanar trajectories. 
Consequently , if the rotational rate of the earth-moon system and the dis- 
tance between the earth center and barycenter are known, then, for a given 
value of V, the associated range of V e may be determined. Through 
the Jacobian relationship, the corresponding value of V at the arrival 
position, point of closest approach to the moon, may be established. 

Since V at the arrival position is perpendicular to the polar position 
vector from the moon center, referred to a nonrotating axis system 

(origin the moon center) may be immediately determined. 

It may be shown that for a given value of C and arrival position 
vector, regardless of the orientation of the position vector, the lunar 
arrival velocity is approximately constant. This assertion must be 
qualified by restricting the altitude to the range considered in this 
study. The justification for the approximation is given in appendix B. 


Determination of Retrograde Velocities 

This section contains a method for determining instantaneous retro- 
grade velocities within certain constrained conditions. For the method 
described herein, the plane of action is defined as the plane containing 
the moon center and the arrival velocity vector. The retrograde impulses 
will be applied in this plane and, therefore, the resulting orbits will 
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be in this plane. Also, since the retrograde impulse is to he initiated 
at the point of closest approach to the moon, then the arrival velocity 
vector will he normal to the position vector. 

The method for determining an orhit with a required pericynthion 
radius and orientation is now presented; the two-body relationships are 
developed in appendix A. The following sketch depicts the vehicle s 
arrival at the vicinity of the moon where the plane of arrival is denoted 
as the X m Y m plane: 



where 0 5 l 80 , measured clockwise or counterclockwise from the r vector. 
% the law of cosines 


V R 2 = V A 2 + V Q 2 - 2V q V a cos 7 


(ID 


where 


cos 7 


rpVop 

rV D 


( 12 ) 


and 



( 13 ) 



9 


Substitution of equations ( 12 ) and ( 13 ). into equation (ll) yields 


S 

1 

0 


GMj 2 * - r) 2 2 V A r p V op 


ar 


V A " - 


It is seen from the relationships in appendix A that 


r p (^ - cos 0) 
a = 


2 — - 1 - cos 9 
r 


It also follows from the relationships in appendix A that 

1/2 

( 1 - cos 0 I 

= |GM| 


'op 


m* 


r pl- - cos 0) 


( 14 ) 


( 15 ) 


(16) 


Substitution of equation ( 15 ) into equation (13) yields 



Now, let 


HI 

,?|x 

0 

(18a) 

GM m = 

A 2 

(18b) 

0 r p 



- COS 


(18c) 


1 - i cos 9 
0 
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- 1) = K 


( l 8 d) 


Substitution of identities (l 8 a) to (l 8 d) into equation ( 17 ) yields 
the following equation: 


v 0 2 = A 2 B + K 


Equation (ll) may now be written as 


V R 2 = V A 2 + A 2 B + K - 2V A 0AB 1/2 (20) 

Equation (20) is the general equation which determines the retrograde 
velocity for the required conditions. 

Transforming equation ( 12 ) yields 


(19) 

S 

1 

0 



By determining 7 then £ can be calculated from the law of 
sines. Thus, 




arc 



(22) 


Equations (20) and ( 22 ) determine the required retrograde velocity 
vector. However, it should be noted that the resultant trajectory could 
represent any type of conic orbit depending on the arrival conditions 
and the required characteristics of the resulting trajectory. It is 
assumed that an elliptical orbit is desired; however, it does not nec- 
essarily follow that the retrograde impulse wilL yield an elliptical 
orbit. The classification of the conic orbit cm be found by determining 
the eccentricity where for an elliptical orbit 0 < e < 1 , for a parabolic 
orbit € = 1 , and for a hyperbolic orbit e > 1 . 

Although it is not envisaged that it will ever be desirable to 
arrive at a certain pericynthion position on a hyperbolic trajectory, 



it is possible to determine the necessary retrograde impulse vector 
which will accomplish this. If equation (20) is modified to 


V R 2 = V A 2 + A 2 B + K + 2V a 0AB 1 / 2 (23) 

then solution of this equation ( 23 ) will yield the necessary velocity. 

The orbit velocity vector in this case will be opposite in direction 
to that which would be determined if equation (20) were utilized. 

In order to obtain an insight into the quantitative values of 
retrograde velocity, refer to figure 4 and table 1. Figure 4 shows the 
retrograde velocities necessary for various combinations of arrival and 
required conditions and the orbit reference table defines the orbits 
established as a result of these retrograde impulses. As an example, 
suppose it is desired to establish an orbit with a pericynthion altitude 
of 100 nautical miles with the pericynthion radius oriented 150 ° from 
the insertion radius. Assume that the vehicle arrives at an altitude 
of 1,000 nautical miles with a velocity of 6,500 ft/sec. Reference to 
figure 4 shows that this will require a retrograde impulse of 
3,360 ft/sec. Reference to table 1 shows that the established orbit 
will have an apocynthion altitude of 1 , 128.6 nautical miles, an apocyn- 
thion velocity of 3 > 031-5 ft/sec, a pericynthion velocity of 
6 , 036.9 ft/sec, and the time required for one complete orbit will be 
3.8 hours. 

It has been determined that the preceding equations will yield the 
retrograde impulse for any specified values of pericynthion radius and 
orientation. However, it is possible in many cases to reach these peri- 
cynthion altitudes with smaller retrograde velocities if the restric- 
tions on the orientation of pericynthion are relaxed. The determination 
of minimum retrograde velocity for a given scalar value of pericynthion 
altitude is given in appendix C. 

A comparison between minimum retrograde velocities and retrograde 
velocities for constrained conditions is shown in figure 5- For the 
construction of this figure, the vehicle is assumed to arrive at an 
altitude of 1,000 nautical miles with an arrival velocity of 7>000 ft/sec. 
Figure 5( a ) shows r p/ r > where r is the arrival radius, plotted against 

0 which will yield the minimum retrograde velocity. Figure 5(b) shows 
r p/ r plotted against the required minimum retrograde velocity and the 

comparison (dashed) curve shows the retrograde velocities which yield 
a pericynthion orientation of 150°. Figure 5(c) shows r p^/r plotted 

against the resultant orbit apocynthion radius for the case where the 
retrograde impulse is a minimum and for the comparison case. It can be 
seen that over a large range of rp/r, the minimum retrograde impulse 
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will offer some fuel savings without influencing the resultant orbits 
significantly. In the range of T ^j T where the fuel savings are con- 
siderable , however , the resultant orbits have the following undesirable 
features: large apocynthion radii, large pericynthion velocities, and 

large orbit periods. 

In order to improve the quantitative feel for retrograde velocities 
as a function of free-coast earth-moon trajectories, retrograde velocities 
for 0 = l80° are plotted against arrival velocity and altitude for a 
fixed earth-moon distance and earth departure position vector in fig- 
ure 6; in sequence (figs. 6(a) to 6 (d)) the value of the departure veloc- 
ity V is varied. On the assumption that the departure angle y is 
between -2° and 25° the normally referred to earth departure velocity 

V e may be quickly determined within the limits of ft/sec and in 
some cases direct reference may be made to figire 3. 


Fuel Requirements 

The fuel required for instantaneous retrograde impulse for estab- 
lishing lunar orbits can be calculated from ths equation 


5 = 1 - 



(24) 


where 6 is the ratio of the fuel mass at commencement of burning to 
the gross mass. Shown in figure 7 are the ratios of fuel mass to gross 
mass necessary for establishing the orbits shown in figure 6. The fuel- 
consumption values as shown are considered to be absolute minimums since 
the impulses are assumed to be instantaneous. Fuel-consumption values 
for corresponding finite impulses, however, will differ little from 
the given values. In all cases where thrust i,5 acting against a 
resolved weight component, there is a loss in efficiency, and this to 
a small degree would be the case with corresponding impulses of finite 
duration. 

As an example of the use of this graph, consider a lunar vehicle 
with an earth-surface weight of 10,000 pounds containing a fuel of 
250 seconds specific impulse. From the graph it will be noted that, 
within the conditions considered, the earth -surface weight of fuel 
required for inserting the vehicle at a lunar altitude of 5,000 nautical 
miles into an orbit with a pericynthion altitude of 100 nautical miles 
is 3 >200 pounds. 


CO rH O 



O H 02 
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In order to obtain the total fuel -requirement for orbit entry and 
exit, the following considerations are taken into account. The retro- 
grade velocity required to insert into the orbit is assumed to equal 
the posigrade velocity for exit. This assumption is made since the 
velocity requirements relative to the center of the moon of the major 
earth-moon trajectory will not have substantially changed after a 
restricted number of lunar orbits. The total fuel requirement is found 
from the equation 


or 



(25) 


6 = 25 - 5 2 


(26) 


where 6 is the ratio of total fuel mass to gross mass required for 
orbit entry and exit. Shown in figure 8 are the ratios of total fuel 
mass to gross mass necessary for orbit entry and exit for the orbits 
shown in figure 6. As an example of the use of figure 8, consider a 
lunar vehicle with an earth-surface weight of 10,000 pounds containing 
a fuel of 250 seconds specific impulse. From figure 8(c) it can be 
deduced that the earth-surface weight of fuel required for both entry 
and exit, at an arrival altitude of 5 >000 nautical miles for an orbit 
with a pericynthion altitude of 100 nautical miles, is approximately 
5,4-00 pounds. 


Orbit Stability 

The simplest and possibly the best method of considering the sta- 
^tlity the vehicle's orbit about the moon is to consider the vehicle 
as a satellite of the moon where the following effects are considered 
to be the major perturbation effects: 

(1) the earth's gravity field 

(2) the sun's gravity field 

(3) the moon's potential distribution 

(4) the lunar librations 

By stability it is implied that successive orbits have repetitive 
characteristics. 
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A brief analysis of the individual gravitational effects of the 
earth, sun, and moon show the possibility of a highly stable orbit. 

The order of magnitude of the gravitational effects of the earth, sun, 
and moon on the vehicle are shown in the following table: 



Gravitational effects, expressed in ft/sec 2 , for - 


Zero lunar altitude 

Lunar altitude of 5,000 nautical miles 

Moon 

5.51 

0.13 

Sun 

.02 

.02 

Earth 

.01 

.01 


If for a short duration there is insignificant difference between 
the effects on the moon and the effects on the vehicle due to the 
gravity fields of the earth and sun, then the moon and vehicle will 
tend to behave as a two-body system with a resultant stable orbit. 

Should the apocynthion altitude of the vehicle's orbit never be greater 
than 5,000 nautical miles, then the scalar accelerations of the vehicle 
and moon due to the sun never have a greater ratio than 1.000147 or the 
inverse. The directional difference in the acceleration vectors is 
negligible since the sun is approximately 80,764,000 nautical miles 
away. The corresponding ratio due to the earth is never greater than 
1.0647 or the inverse and the directional difference is small. Viewed 
in this manner, the sun's gravitational effect is very small. The moon 
and vehicle acceleration vectorial difference due to the earth appears 
to be more significant, but due to the oscillatory nature of the dif- 
ference (since it is periodic with the vehicle's orbit), the effect is 
small over a restricted number of orbits . 

Other perturbation effects of interest are the librations of the 
moon about its center of gravity. The apparent librations as viewed 
from earth are of no concern in this study. There is a small real 
libration in longitude due to the eccentricity of the moon's orbit 
about the barycenter but the period is a month, and hence this libra- 
tion does not present a problem. If any high-frequency librations 
exist, they are thought to be insignificantly small. 

The main changes in orbit characteristics with time are anticipated 
to be changes in the inclination of the orbit and regression of the nodes 
relative to the lunar equator. 

Several stability checks have been conducted by a comprehensive 
simulation incorporated in a digital mathematical model which includes 
earth, sun, and lunar potential distribution effects. 




This model is particularly attractive in that the origin of integra- 
tion is the moon for the conditions considered and the round-off errors 
are those involved in the integration of perturbations from reference 
ellipses. The checks added confidence to the preceding discussion of 
stability. 


CONCLUSIONS 


In a preliminary study of retrograde velocities required for 
insertion into lunar low-altitude orbits, the following conclusions 
were reached: 

1. At low lunar altitudes, say less than 5,000 nautical miles, the 
lunar vehicle in free-coast conditions will essentially behave as a 
satellite of the moon and, hence, the trajectory will be conical with 
the center of the moon mass as a focal point. Depending on the velocity 
imparted by the retrograde impulse, the trajectory will be elliptic, 
parabolic, or hyperbolic. Elliptical trajectories are of interest in 
that closed orbits about the moon are required for survey purposes. 

2. A simple thrusting logic may be introduced for obtaining a 
required pericynthion position, the only requirement being that the 
insertion velocity is that which would be yielded by the classical two- 
body solution for the required conditions. 

3. The lunar arrival velocity increases with increase in earth 
insertion velocity and this entails heavier fuel expenditure for inser- 
tion into a local lunar orbit; however, there is the possibility that 
the trajectory associated with a higher earth insertion velocity will 
require less fuel expenditure for guidance before the arrival phase. 

b. For a given arrival altitude, the nearer the required orbit is 
to a circular orbit, the lower the required retrograde velocity. Further- 
more for a given pericynthion altitude, the nearer the associated orbit 
is to circular, the smaller the pericynthion velocity. 

5- For an elliptical orbit with pericynthion altitude and orienta- 
tion stipulated, there is a unique retrograde velocity. If the restric- 
tion on the orientation of the pericynthion is relaxed, then in many 
cases it is possible to determine a smaller retrograde impulse which 
will yield the desired pericynthion altitude. 
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6. There are cases when the resultant perlcynthion of a minimum 
retrograde impulse does not occur directly opposite the point of appli- 
cation. Utilization of a minimum retrograde impulse in these circum- 
stances may produce an orbit with an excessive apocynthion altitude, 
perlcynthion velocity, and period. It is possible to minimize these 
adverse features, at the expense of fuel, by orienting the perlcynthion 
away from the point of minimum retrograde impulse to a point where 
the orbit characteristics become more desirable. 


Space Task Group, 

National Aeronautics and Space Administration, 
Langley Field, Va. , June 12 , 1961. 



APPENDIX A 


TWO -BODY EQUATIONS 


In order to maintain continuity In the body of the text, transposes 
of well-known two-body relationships are immediately used. This appendix 
is included to indicate the derivation of these relationships. 

The differential equations of motion are 


and 


2 

dt 2 



®yn 

r 2 


m 

r 



= 0 


(Al) 


(A2) 


From these equations and a knowledge of conic geometry, the following 
equations can be evolved. The general conic expression for the semilatus 
rectum is 

p = r( 1 + e cos 9) (A3) 

and when r = r p , then 0=0 which yields 

p = r p (l + e) (Ak) 

The expression for the angle y which is the angle between the 
velocity vector V Q and the normal to the vector r is from integra- 
tion of equation (A2) 


cos y = 


r p^op 

rV 

v o 


(A5) 


The expression for the semimajor axis which is evolved from equa- 
tion (A3) is 
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1 - 6 " 


Substitution of the values of p and € from equations (A3) and (a 4) 
yield 


r P r “ COS 9 


2 — - 1 - cos 0 
r 


For a hyperbola since € > 1 it can be seen from equation (A6) that the 
numerical value of a will be negative. The orbit velocity is 


V 0 2 « GMJ§ - I 


Substituting equation (AT) into equation (AS) yields 


_ r / 1 - cos 


P 1 . JL cos 0 




The pericynthion velocity is 


Vcp 2 * 


Substitution for a from equation (A7) in equation (A10) yields 


= GMm 


1 - cos 8 

rv/— - ccs 0 
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APPENDIX B 

ARRIVAL VELOCITY APPROXIMATIONS 


S 

1 

0 


This appendix is presented to show that for a given value of C 
the arrival velocity is almost entirely dependent on lunar altitude, 
within the limits of this study. This is succeeded by a derivation of 
an approximation of the variation of arrival velocity with lunar altitude. 


Arrival Velocity at Constant Lunar Altitude 
The Jacobian velocity relationship is as follows: 


„2 2 2 . 2(l - n)GM 2 hGM 

V = “ ^bp + = - + — ~ — 

-^ep r 


(Bl) 


For convenience, equation (Bl) is rewritten so that the terms independent 
of the orientation of r appear on the left-hand side as 


y2 _ 2 |_iGM 
r 



+ 


2(1 - [jl)GM 
r ep 


(B2) 


Upon utilizing the law of cosines, the following equation is evolved: 


^bp 2 


2(1 - q)GM 


-ep 


= 0)2 


(R m 2 + r2) + 2R m ( X - R^ 


+ 


2(1 - n)GM 






1/2 


(B3) 



20 


Since 






2 + r 2 ) »> 


2Rm(x - Rm) 


and. 


\ »> r; X - R m < r 


+ 


+ r 


»> 2 (^m + ?e)(x - Rm) 




Therefore, 


^ 2 + 2(1 - M-)GM 


+ r 


ep 


CD 


(Rm 2 + r 2 ) 


2(1 - pi)GM 


(Rm + Re) 2 + ^ 


1/2 


This allows equation (Bl) to be modified to a close approximation for 
the range of altitudes studied as follows: 


V 2 « CD 2 (R m 2 + r 2 ) + 2 U - K) GM ... + - c (b4) 

(Rm + Re) + ^ 


s 

1 

0 


Consequently, it can be seen from approximation (B4) that if r is 
constant, then the orientation of r will non affect the value of V, 

— > — > f — > — A — ) ^ 

V m = Y + 1 a) X r I where a) is a constant vec:or and r has a constant 

scalar value. If trajectories in the earth-moon plane are being con- 
sidered, then V m = V ± cor where cor is a constant. 


Variation of Arrival Velocity With junar Altitude 

The trend of the arrival velocity with increase of altitude above 
the lunar surface may be more clearly understood by the following con- 
siderations. Substitution of equation (B3) iato equation (Bl) yields 
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S 

1 

0 


V 2 = o£|R m 2 + r 2 + 2R m (X - R, 


^m 


2(1 -,u)GM 


(Re + Rm) 2 + 2 (R e + R m )(x - R m ) 


+ C 


1/2 


(B5) 


Differentiation of this equation with respect to r yields 


SV 2 

Sr 


= 2u? 


r - 


2(1 - n)(GM)r 


(R e + ^m^ + r2 + 2 (r 6 + Rm^(x - R^ 


2|_lGM 


5/2 r2 


(b6) 


however > 

r 




l 2a£z 


2(1 - n)(GM)r 


«< 


2|_iGM 


V. 


^R e + R m ^ + -i- 2 ^R e + 


1 5/2 


Therefore, it is possible to reduce equation (B6) to the approximation 


5 v 2 _ -3.5 x io l4 

Sr p2 

This expression gives a good approximation for the trend of the arrival 
velocity with altitude. Although this approximation applies for a 
velocity V, which is referred to the rotating axis system with origin 
at the barycenter , it is now shown that with little loss of accuracy V 
may be regarded as V A • Given that 


ft/ 


sec 1 


(B7) 


(B8) 


V A 2 = V 2 ± 2curV + o^r 2 
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then. 


dV 


A 

dr 


!r + * ^ 2mrV) 


(B9) 


which yields 


dV 2 cur Sv 2 o 

~ sr 1 t sr + 2<A ' 1 ^ 


(BIO) 


S 

1 

0 


and upon collecting terms 


3V/i 


- *!/i ± * 


Sr Sr y 


2co(<jur ± V) 


( Bll) 


^at nSSrtin ^ a relevant ran 8 e of values into ihe equation it is found 


cor 

T <«! 


and 


i 2 

| 2CD( 6DT ± V)|<« ^ L 
Sr 


Therefore, the following approximation is acceptable: 


^A 2 ^ ^ -J.5 x lO 1 ^ 

Sr Sr o 


ft/sec^ 


(B12) 
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APPENDIX C 

DERIVATION OF THE MINIMUM RETROGRADE VELOCITY TO 
OBTAIN A GIVEN PERICYNTHION SCALAR VALUE 

Upon referring to the following equation 

V R 2 = V A 2 + A 2 B + K - PV^AB 1 / 2 (Cl) 

and the identities (l8a), (l8t) , and (l8d), it is seen that the quan- 
tities A, K, and 0 are independent of 0; therefore, differentia- 
tion of equation (Cl) with respect to 9 yields 



and differentiation of identity (l8c) yields 



(C2) 


(C3) 


A necessary condition for minimum retrograde velocity is that 
and it is known that A / 0; therefore, 


3Vr 2 

89 


0 , 




2k 


Suppose, tentatively, that for equation C4) to hold, 


A - V A 0B _1 / 2 = 0 


(C5) 


therefore ; 


B' 


-1 


(V A ! 0) £ 


Substitution of equation (C5a) into identity (l8c) yields 


(C5a) 


S 

1 

0 


cos 9 = 




i / a y 
0 " 1va0) 


(c6) 


Equation (C4) has two boundary conditions 9=0° and 9 = l80°. 
When r = r p; then 0 = 0 by definition of tie two-body orbit and this 

is the unique value of 9 to be considered. By reference to equa- 
tion (C4), it is seen that 9 = l80° is an alternative solution to 


cfV 


R 


1 - “\sin 9 

0/ 


89 


= 0 


1 - 


cos 9 


(C7) 


Two conditions are now stipulated. For condition 1 where 9 = l80° 
and ll— = o let Vr = Vpq. For condition 5 where l80 = j 9 | ^ 0 


89 


and 




R 


8b 


0 let = V ^2 ■ 


The following analysis is to prove that when V^? exists , it is a 

2 2 

minimum. Tentatively assume that ^ Vp£ 


in which case the 
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S 

1 

0 


following equation is evolved from equation (Cl): 


A 2 B x - 2V a 0AB 1 1 / 2 > A 2 B2 - 2V a 0AB2 1//2 


(c8) 


It is convenient to add 2V A 0ABpl/ 2 - A 2 B 2 to both sides. This 
allows inequality (c8) to be rewritten as 


(bi L/2 + ^ 2 1 / 2 )( b 1 1/2 - B2 1/2 ) £ 2V a 0A^! 1/2 - Bg 1 / 2 ) 


( C8a) 


^ „ 1 - cos 9 

By definition, B = ; therefore, 

1 - i cos 9 

0 


Bl 


1 / 2/2 


1/2 


n + i 


0 


(C9) 


For condition 2 equation (C6) yields 


1 - 


-1 < 


/ A \ 

W 


1 


A 

yl? 


(CIO) 


By introducing the value of 


(CIO) simplifies to 


\ A0 , 


from equation (C5a), inequality 


B2 ? 


> 


1 + i 


(Clla) 
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and 


Ba 1/2 > ' 2 


, 1/2 


1 + 


0 / 


( Cllb) 


From equation (C9) and inequality (Cllb) it follows that 


- B2 1//2 < 0 


(C12) 


S 

1 

0 


Inequality (c8a), which is the condition for V R1 = V R2 which 
in turn implies that is a fldnimuin, is now divided by 

,1/2 _ t^1/2)a2 


a negative quantity which reverses the inequality, 

and hence 


Bj 1 ' 2 + B2 1 / 2 < ^ 


(ci3) 


From equation (C5a), both sides of which are positive, it is seen 


that 




(C14) 


Therefore, substitution of equation (Cl4) into inequality (C13) yields 


Bl 


!/ 2 + B2 1 / 2 S 2B2 1/2 


(C15) 


l/2 

Subtraction of 2 B 2 ' from both sides of inequality (C15) yields 


B^ 2 - Ba 1 / 2 S 0 


(C16) 
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but from inequality (C 12 ), it is known that this is true: therefore 

uhPn 4- ^ j 4- -r „ _ • _ 7 J 


it is a minimum. It must be noted that V D o is 

Kc 


when Vj^2 exists, 

minimum for a conic trajectory. For the trajectory to be elliptical, 
two-body orbital equations yield equation (15) and for elliptical con- 
ditions to exist 


S 

1 

0 


or 



cos 0 > 0 



1 > cos 0 


(C17) 


Hence, when Vro exists and 2 -2 _ 1 > cos 9 thpn v i. +k# . 

^ r x,nen is the minimum 

retrograde velocity for an elliptical orbit. 

In summation, if V Rg does exist, it is a minimum and, of necessity, 

. , 1 " (w 


80 


= 0 and -1 < 


< 


1. 


i- U. 

* V 


tion (C 4 ). 


When 0 - 1 , 0-0 is the unique value which satisfies equa- 


If Vr 2 does not exist and 001 , then V R1 is a minimum. If 


V R1 is a minimum, 0 = l8 0 ° and, of necessity, ** - = 0 . also implied 


80 


is that 2 — - 1 > 
r 


cos 0, which is the condition for ellipticity. 


The method for determining the minimum retrograde 
given set of conditions is outlined as follows: 


velocity 


for a 



28 


Step 1: 

Determine from the following equation if does exist: 


(Cl8) 

S 
1 
0 

If there is a solution to equation (Cl8), then V R2 does exist and the 

determined value of 0 will yield the minimum retrograde impulse. 

Should the equation have no legitimate solution, then Vp2 &oes not 

exist and, therefore, 0 = l80°, which is the value of 0 for Vp>p, 
will yield the minimum retrograde impulse. 

Step 2: 

The eccentricity of the orbit should now be determined. By using 
the value of 0 as found from equation (Cl8), the eccentricity can be 
determined from 

e = j - 1 (C19) 

Step 3: 

The minimum retrograde velocity for an elliptical orbit can be 
ascertained by using the determined value of 0 In equation (20) which 
is 

V R 2 = V A 2 + A 2 B + K - 2V a 0\B 1 / 2 

If the resultant orbit is not elliptical, the minimum retrograde velocity- 
can be determined from equation (23) which is 

Vr 2 = V A 2 + A 2 B + K + 2V A j2iAB 1 / 2 
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TABLE 1 


ORBIT REFERENCE 



9 = 120° 

9 = 150° 

9 - l8o° 

T 

| Arrival altitude, 1,000 nautical miles 

L _____ _ 

: For hp =» 50 nautical miles, 

ha, nautical miles . . . . 

Voa, ft/sec 

j V G p, ft/sec 

Period, hr 

1,914.0 

2,263.8 

6,535.5 

5.3 

1,143.6 

2.963.8 

6.244.9 
3-8 

1,000.0 

3,146.1 

6,171.7 

3-5 

For hp = 100 nautical miles, 
ha, nautical miles .... 

Voa j ft/sec 

V op , ft/sec 

Period, hr 

1,788.0 

2,397.1 

6,296.8 

5.1 

1,128.6 

3,031.5 

6,036.9 

3.8 

1,000.0 

3,197.5 

5,970.9 

3.6 

For hp = 200 nautical miles, 
ha, nautical miles ... . 

Voa, ft/sec 

V 0 p, ft/sec 

1 Period, hr 

i . . .. 

1,593.0 

2,638.5 

5,868.9 

4.9 

1,102.9 

3.156.7 

5.661.8 
4.0 

1,000.0 

3,293.2 

5,608.7 

3.8 

Arrival altitude, 3,000 nautical milei 

For hp = 50 nautical miles, 
ha, nautical miles . . . . 

V^, ft/sec 

V 0 p, ft/sec 

Period, hr 

a 844,339-7 

8.8 

7,578.9 

17,330.2 

4,C f4.3 
1,366.0 
6,551.0 
L0.3 

3,000.0 

1,701.0 

6,780.4 

7.7 

For hp - 100 nautical miles, 
ha, nautical miles .... 

Voa, ft/sec 

Vop , ft/sec 

Period, hr 

*56,562.3 

132-3 

7,333.1 

315.3 

3,588.0 

1,416.6 

6,724.4 

10.2 

3,000.0 

1,734.6 

6,582.3 

7.8 

For hp - 200 nautical miles, 
ha, nautical miles .... 

Voa, ft/sec 

V D p, ft/ sec 

Period, hr . 

*20,965.9 

357-9 

6,889.2 

77.9 

3,e+5-l 

1,5 L0.8 
6,350.8 
L0.1 

3,000.0 

1,798.3 

6,223.8 

8.0 

Arrival altitude, 5,000 nautical miles 

For hp * 50 nautical miles, 
ha, nautical miles .... 

Voa, ft/sec 

Vop, ft/sec 

Period, hr 

Hyperbolic 

Trajectory 

8,^55. 5 
167-9 
7,2 39*1 
23-1 

5,000.0 

1,168.2 

7,021.7 

12.8 

For hp = 100 nautical miles, 
h a , nautical ml lee .... 

Voa, ft/sec 

V 0 p, ft/sec 

Period, hr 

Hyperbolic 

Trajectory 

8, c 44.8 

1 09.3 

7,C05-4 

22.3 

5,000.0 

1,193.0 

6,826.5 

13.0 

For hp = 200 nautical miles, 
ha, nautical miles .... 

Voa, ft/sec 

V 0 p, ft/sec 

Period, hr 

Hyperbolic 

Trajectory 

7,! 78.9 
836.6 
6 ,( 36.5 
21.1 

5,000.0 

1,240.4 

6,472.8 

13-2 


a In the solar system the assumption of elliptical characteristics in these 
cases is not valid. 
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Application of posigrade impulse 



Figure 1.- Diagram of envisaged manned lunar mission. 





ft/sec 







Figure 6.~ Retrograde velocity for insertion into orbit for 9 = l8o°. Earth departure 
conditions; altitude, 100 nautical miles; moon lead angle, 90 ° > earth-moon distance, 
192,358 nautical miles. 
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(b) V = 35,796 ft/ sec . 


Figure 6.- Continued 




(c) v. = 35,885 ft/sec. 
Figure 6.- Continued. 
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(d) V = 35, 97**- ft/sec. 


Figure 6.- Concluded. 
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Figure 7.- Fuel consumption for insertion into orbit. 
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